Accurate Test of Chiral Dynamics in the jp — > ir°p Reaction 
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A precision measurement of the photon asymmetry E and differential cross sections da/dQ, for 
the ■yp — > 7r°p reaction in the near-threshold region has been performed with a tagged photon beam 
and almost 4-7T detector at the Mainz Microtron. The Glasgow-Mainz photon tagging facility along 
with the Crystal Ball/TAPS multi-photon detector system and a cryogenic liquid hydrogen target 
were used. These data allowed for a precise determination of the energy dependence of the real 
parts of the S- and all three P-wave amplitudes for the first time and provide the most stringent 
test to date of the predictions of Chiral Perturbation Theory and its energy region of convergence. 
The upper limit of agreement is between 165 and 175 MeV incident photon lab energy, ~ 25 MeV 
above threshold. 



Low-energy pion photoproduction experiments are of 
special interest because the pion, the lightest hadron, is a 
Nambu-Goldstone boson that by its existence represents 
a clear signature of spontaneous chiral symmetry break- 
ing in QCD ill]. The dynamic consequences are that the 
production and elastic scattering of neutral pions at low 
energies are weak in the S-w&ve and strong in the P- 
wave [U-1H i as is seen in the 7^V — > nN reaction 0, 0] . In 
neutral-pion photoproduction the S'-wave amplitudes are 
small since they vanish in the chiral limit (m„, — > 0); 
their small, but non- vanishing values are measures of ex- 
plicit chiral symmetry breaking. In addition, they are 
isospin violating 0, H, 0| since m u ^ md [a 0- The 
magnitudes of low-energy scattering and production ex- 
periments are predicted by Chiral Perturbation Theory 



(ChPT), an effective field theory of QCD based on spon- 
taneous chiral symmetry breaking [lH5( . Our efforts have 
been focused on accurate measurements of low-energy 
7./V — > -kN reactions, including the sensitive spin observ- 
ables, to perform tests of these predictions. As has been 
stressed [10|, any serious discrepancy between these cal- 
culations and experiment must be carefully examined as 
a potential violation of QCD. 

In 2008 we conducted an investigation of the jp — > TT°p 
reaction with the twin goals of obtaining 1) the energy 
dependence of the photon asymmetry E for the first time, 
and 2) the most accurate measurement to date of the dif- 
ferential cross section from threshold through the A re- 
gion. The energy dependence of E, in combination with 
the cross-section data, allowed us to extract the real parts 
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of all P-wave as well as the S'-wave multipoles as a func- 
tion of photon energy in the threshold region. These 
data in turn also allowed the first test of how well ChPT 
calculations agree with the data as a function of photon 
energy above threshold. There exists one previous mea- 
surement of the photon asymmetry |ll| , but due to poor 
statistics resulting from small cross sections and limited 
detector acceptance (fts 10% for tt° detection) £ was in- 
tegrated over the entire incident photon energy range, 
leading to data at only the bremsstrahlung-weighted en- 
ergy of E 1 = 159.5 MeV. Moreover, the contribution 
to the asymmetry from the target walls — significant in 
the threshold region — was not properly taken into ac- 
count [HI ■ With the present set-up the azimuthal accep- 
tance was vastly superior and symmetric, the degree of 
linear polarization was higher, a rigorous empty-target 
subtraction has been done, and as a result both the sta- 
tistical and systematic uncertainties are much smaller for 
£ as well as the cross sections. The most accurate pre- 
vious measurement along with a list of earlier efforts can 
be found in Ref. 



The experiment that is the focus of this letter was 
conducted at the Mainz Microtron MAMI [H, [l4| 
where linearly polarized photons, produced via coher- 
ent bremsstrahlung in a 100-/im-thick diamond radia- 
tor [H EH, were sent through a 4-mm-diameter Pb 
collimator and impinged on a 10-cm-long liquid hydro- 
gen (LH2 ) target located in the center of the Crystal 
Ball 1171 . The TAPS spectrometer served as a forward 
wall [18( , and the LH 2 target was surrounded by a particle 
identification detector [19| , used to differentiate between 
charged and uncharged particles. The incident photons 
were tagged up to an energy of 800 MeV using the post- 
bremsstrahlung electrons detected by the Glasgow-Mainz 
tagger [20]. For the electron beam of 855 MeV used in 
this experiment, the tagger channels had a width of about 
2.4 MeV in the tt° threshold region. The diamond radia- 
tor was positioned relative to the electron beam such that 
the photons produced had a polarization in the range 50- 
70% between ir° threshold and 
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200 MeV 

Neutral pions produced in the LH 2 target were iden- 
tified in the detector system using their 27 decay 
and a kinematic-fitting technique described in detail in 
Ref. [2l|. Both two- and three-cluster events that satis- 
fied the hypothesis of the reaction jp ir°p — > 77P with 
a probability of more than 2% were accepted as candi- 
dates for this reaction. Background contamination of the 
event candidates was found to be from two sources: inter- 
actions of the bremsstrahlung photons in the target ma- 
terial different from liquid hydrogen, and accidental co- 
incidences between the tagger hits and the trigger based 
on the detector signals. The background was subtracted 
from the signal directly by using two different data sam- 
ples, the first of which included only events with acci- 
dental coincidences, and the second was taken with an 
empty target cell. 



Acceptance of the detector system was determined by 
Monte Carlo simulation of jp — > ir°p using an isotropic 
angular distribution. All events were propagated through 
a GEANT3.21 simulation of the experimental set-up, 
folded with resolutions of the detectors and conditions 
of the trigger. Close to the reaction threshold, the 
production-angle acceptance was found to be almost uni- 
form with a detection efficiency about 80%. 

The systematic uncertainties in the absolute numbers 
of the differential cross sections for the reaction jp — > tt°p 
obtained in the analysis of the data set were estimated to 
be not larger than 4%. Such a magnitude of the system- 
atic uncertainty is mostly determined by the calculation 
of the photon-beam flux, the experimental detection ef- 
ficiency, and the number of protons in the LH 2 target. 
The systematic uncertainties in the numbers for the pho- 
ton asymmetry are on the level of 5%, where this value 
comes mostly from the uncertainty in the determination 
of the degree of the linear polarization of the incoming 
photons. 

Results for the differential cross section and the pho- 
ton asymmetry are presented in Fig. [1] as a function of 
pion center-of-mass (CM) production angle 9 at E 1 — 
163.4 ±1.2 MeV, and as a function of incident photon 
energy at 8 ~ 90 ± 3°. Also shown are three theoreti- 
cal curves: 1) Heavy Baryon Chiral Perturbation Theory 
(HBChPT) calculations to 0{q A ) [H with the five em- 
pirical low-energy constants brought up to date by fitting 
these data [23]; 2) relativistic ChPT calculations (also to 
0(q 4 )) which as well have five low-energy constants fit 



to these data [24j, l25j; and 3) an empirical fit with error 
bands calculated using the formalism from Refs. 26, 27j j . 

Due to the high quality of the present data, it is, for 
the first time, possible to determine the energy range for 
which ChPT agrees with the data. The values of the low- 
energy constants were obtained from fits to the data in 
the range from 150 MeV to a variable P™ ax . This was 
done 1231 using the 0(q A ) formulas of heavy baryon calcu- 



lations [22] and also for the relativistic theory [24ll25l.l28|. 
Figure 0Je) displays the reduced \ 2 f° r the empirical fit 
and both ChPT calculations. For P™ ax up to ~ 167 MeV 
the ChPT calculations are consistent with the empiri- 
cal fit, but above this energy the relativistic calculation 
starts to deviate from the data whereas the heavy baryon 
calculation begins to deviate at ~ 172 MeV. This is 
interesting since the relative contributions of the terms 
containing the low-energy constants of the relativistic cal- 
culation are significantly smaller than those for the heavy 
baryon version, suggesting a better convergence for the 
relativistic ChPT method. 

The next step in the interpretation of the data is to 
extract the multipole amplitudes and compare them to 
the theoretical calculations. To set the notation, the dif- 
ferential cross sections can be expressed in terms of the 
S- and P-wave multipoles (E + , Pi, P 2 , P3) and can be 
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FIG. 1. (color online) Differential cross sections in /ib/sr (a) 
and photon asymmetries (b) for 7r°-production as a function 
of pion CM production angle 6 for an incident photon energy 
of 163.4 ± 1.2 MeV. Energy dependence of the differential 
cross sections (c) and photon asymmetries (d) at 6 ~ 90 ± 3° . 
Errors shown are statistical only, without the systematic un- 
certainty of 4% for da/dQ. and 5% for E. The theory curves 
are dashed (black) HBChPT (2jJ, dash-dotted (blue) rela- 
tivistic ChPT 25, 28], and (green) an empirical fit with an 
error band, (e) Reduced y 2 as a function of the incident pho- 
ton energy for HBChPT [23j (open black circles), relativistic 
ChPT [25II28I] (open blue triangles) and an empirical fit (solid 
green dots), with lines drawn through the points to guide the 
eye. Note that in (c) and (d) the two points in incident pho- 
ton energy below ir + threshold are included; these two points 
are excluded in the fits shown in (e) due to their large error 
bars. 



written as 



BcosO + Ccos 2 



(1) 



where q and k denote the CM momenta of pion and 
photon, respectively. The coefficients are given by 
A = \E Q+ f+P 2 l with P 2 1 = KW + W), B = 
2Re(E +P?), and C = |Pi| 2 - P 2 %. The measurement of 
the cross sections of earlier experiments [ll[ permitted 
the extraction of E + , Pi and the combination P23 • In 



order to extract the values of ReE + and all three P- 
waves separately from the data, it is necessary to also 
measure the photon asymmetry 



d<7± — dan 
da± + dan 



«/§(«). (2) 



where da± and da^ are the differential cross sections for 
photon polarization perpendicular and parallel to the re- 
action plane with pion and outgoing proton. To reiterate, 
the measurement of the differential cross section and S 
allows for the separation of the four multipoles. It is im- 
portant to note that the determination reported here is 
more accurate than previous ones due to the far smaller 
uncertainties of the cross sections as well as the energy 
dependence of £. Furthermore, we note that the D- waves 
have been neglected in both ((T|) and ^ , but they have re- 
cently been shown to be important in the near-threshold 
region [29|, [3(| . Since there are insufficient data to deter- 
mine the D-wave multipoles empirically, they have been 
taken into account by using their values in the Born ap- 
proximation, which is sufficiently accurate for the present 
analysis. 

The empirical fits to the data start with the following 
ansatz for the multipoles 



E 0+ (W) = E® + £g? ( ^— ^ ) + , (3) 
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where here E 1 and E^ r are in the lab frame 
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P^ (with i 



and 

1,2,3) are constants that 



are fit to the data. Based on unitarity, the cusp 
parameter in has the value = ReE +(jp — > 
7r + n)a ccx (7r + n — > Tr°p) [31]. Using the best experimen- 
tal value of a cox (7r-p -> 7r°n) = -(0.122 ± 0.002)/?^+ 
obtained from the observed width in the Is state of 
pionic hydrogen (32I ] . assuming isospin is a good sym- 
metry, i.e. a ccx (-!T + n -t ir°p) = -a ccx (ir~p -)• 7r°77), 
and the latest measurement for Ep+ (jp ir + n) = 
(28.06 ±0.27 ±0.45) x 10" 3 /m T + 0, we obtain /3 = 
(3.43 ±0.08) x 10^ 3 /m^+, which was employed in the 
empirical fit. The uncertainty introduced by the errors 
in /3 and isospin breaking are smaller than the statisti- 
cal uncertainties of the multipole extraction depicted in 

Fig. m 



1 The empirical values are in units of 10 ^ jm^ 
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The extracted multipoles are displayed in Fig. [2] along 
with the theoretical calculations. The points are single- 
energy fits to the real parts of the S- and P-wave mul- 
tipoles, and the energy-dependent fits from © and (U]) 
are shown with the error band. The imaginary part of 
the S'-wave multipole E + was taken from unitarity ([3]) 
with the value of the cusp parameter explained above, 
the imaginary parts of the P- waves were assumed to be 
negligible, and the P-wave multipoles were calculated in 
the Born approximation. The impact of P-waves in the 
P-wave extraction is negligible [29[ but in the S'-wave it 
can be sizeable. In order to assess the uncertainties in the 
S-wave extraction associated to our P-wave prescription 
we have estimated the uncertainty from the difference be- 
tween the Born terms and DMT dynamical model in (34| . 
This error estimation is depicted in Fig. [5] as a gray area 
at the top of the first plot. 

As was the case for the observables, there is very good 
agreement between the two ChPT calculations and the 
empirical values of the multipoles for energies up to ~ 
160 MeV with the same pattern of deviations above that. 

In conclusion, the combination of the photon asymme- 
try and improved accuracy in the differential cross sec- 
tion has allowed us to extract the real parts of the S-wave 
and all three P-wave multipoles as a function of photon 
energy for the first time. As was the case for the observ- 
ables, the ChPT calculations agree with the multipoles 
up to an energy of ~ 165 to 175 MeV, with the relativis- 
tic calculations deviating at the lower-energy end and the 
HBChPT calculations closer to the upper end. 
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